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ABSTRACT 

Results are presented of a new VLA— ROSAT study that probes the magnetic field strength and 
distribution over a sample of 16 "normal" low redshift {z < 0.1) galaxy clusters. The clusters span two 
orders of magnitude in X-ray luminosity, and were selected to be free of (unusual) strong radio cluster 
halos, and widespread cooling flows. Consistent with these criteria, most clusters show a relaxed X-ray 
morphology and little or no evidence for recent merger activity. 

Analysis of the rotation measure (RM) data shows cluster-generated Faraday RM excess out to 
^0.5 h^l Mpc fr om cluster centers. The results, combined with RM imaging of cluster-embedded 
sources and ROSAT X-ray profiles indicates that the hot intergalactic gas within these "normal" 
clusters is permeated with a high filling factor by magnetic fields at levels of < \B\ >icm = 5—10 

{£/lO kpc)^^/^ /lyg^ /iC, where £ is the field correlation length. These results lead to a global estimate of 
the total magnetic energy in clusters, and give new insight into the ultimate energy origin, which is likely 
gravitational. These results also shed some light on the cluster evolutionary conditions that existed at 
the onset of cooling flows. 

Subject headings: galaxies: clusters: general — magnetic fields — polarization — radio continuum: galaxies — X-rays: 
general 



1. INTRODUCTION 

Determinations of magnetic field strengths in the intra- 
cluster medium over the past decade have revealed fields 
of unanticipated strength in some clusters. Additionally, 
these studies have raised apparent contradictions between 
field strengths (or limits) obtained from different observa- 
tional methods. In order to understand the physical condi- 
tions in the intracluster medium, it is important to under- 
stand the origin(s) of cluster fields, and the astrophysical 
processes that relate them to other energy constituents of 
the cluster gas. 

Estimates of magnetic field strengths in the intracluster 
medium (ICM) can be determined from measurements of 
Faraday rotation through the intracluster gas, and an in- 
dependent measurement of the ICM thermal electron den- 
sity. He- Faraday rotation is given by: 



RM = ^ = 811.9 



n,Bn 



rad 
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AA2 ^ ^0 ^ 

where x is the position angle of the linearly polarized ra- 
diation at wavelength A, is the thermal electron density 
in cm~^, is the line of sight magnetic field strength in 
fiG, and L is the path length in kpc. The thermal elec- 
tron density in a cluster can be determined from X-ray 
surface brightness profiles of the hot (T~ 10^ K), diffuse 
(n 

the cluster potential 

The first study of background RMs over a single clus- 



^75^ cm gas (Bohringer 1995) which fills 



ter (Kim et al. 1990) was a targeted set of deep VLA 
observations of 18 sources close in angular position to the 
Coma cluster (which has no cooling flow, but a strong syn- 
chrotron halo). The result of the study was < \B\ >icm = 

2.5 {£/10 kpc)^^/^ /lyg^ /uG where £ is the B-reversal scale. 
In a subsequent study by Feretti et al. (1995) the discov- 
ery of smaller £ scales down to 1 kpc raised this estimate 

to ~ 7.2 ft,75^/iG for Coma. 

Studies of a large number of galaxy clusters with many 
RM probes per cluster are currently unfeasible given the 
sensitivity limits of available radio telescopes, combined 
with the small angular size of more distant clusters. This 
situation can be circumvented by obtaining RM probes 
through a sample of clusters, each having typically only 
one or two (bright) polarized radio sources. The consen- 
sus of these studies (Lawler & Dennison 1982; Kim et al. 
1991; Goldshmidt & Rephaeli 1993) is that cluster cores 
have a detectable component of RM, and many have fleld 
strengths at the microgauss level. Note that the Faraday 
study by Hennessy, Owen, & Eilek (1989) does not flnd 
evidence for intracluster magnetic fields. This discrepancy 
is, however, likely due to the combination of small statis- 
tics and large impact parameters in their study. 

For a few clusters with extensive cooling fiows, high 
resolution Faraday rotation measure mapping of extended 
sources that are embedded within the cooling fiow zones 
have, in combination with X-ray data, produced magnetic 
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field strength estimates of 10-40 fiG, ordered on scales 
varying from 100 - 0.5 kpc (see Taylor, Allen, & Fabian 
1999, and references therein). 

Lower limits to ICM magnetic field strengths in the 
range 0.1 1 fiG level have been suggested from recent de- 
tections of both excess (over thermal) extreme ultraviolet 
(EUV) and hard X-ray (HEX) emission in some clusters. 
These field values would seem to prima facie contradict 
the much higher values above if the EUV excess emission 
is interpreted as inverse Compton (IC) scattering of ^ 100 
MeV electrons (Rephaeh, Ulmer, & Gruber 1994). The 
EUV detections (Bowyer, Berghofer, & Korpela 1999) in 
particular apply to only a few clusters having widespread 
extended synchrotron emission and/or they occur near to 
the cooling flow zone (as in M87). Spatial differentiation of 
high and low magnetic field regions in the ICM avoids the 
apparent contradiction by allowing the synchrotron and 
EUV IC emission to originate in low field regions, while 
high field regions, where the synchrotron energy loss time 
is short, would provide the major contribution to the Fara- 
day rotation measures (Enl31in, Lieu, & Biermann 1999). 
Further, EUV emission also appears not to be cluster wide, 
but concentrated to central sub-regions, and hence does 
not spatially correspond to the wider regions probed by 
the X-ray emission and Faraday RM measurements. The 
HEX excess in clusters can be plausibly understood as 
bremsstrahlung from a (probably shock heated) popula- 
tion of suprathcrmal electrons (Enfilin, Lieu, & Biermann 
1999; Dogiel 1999; Sarazin & Kempner 2000). Thus, 
on balance, it appears that spatial differentiation of field 
regions and/or emission mechanisms other than IC scat- 
tering of CMB photons are more plausible where apparent 
contradictions in magnetic field estimates arise. 

This paper concentrates on what we shall term "normal" 
Abell clusters i.e. those which have neither widespread 
cooling fiows nor strong synchrotron halos. The obser- 
vations are aimed at estimating the strength and spatial 
extent of intracluster medium magnetic fields for a rela- 
tively homogeneous sample of 16 Abell clusters. Through- 
out this Letter, we adopt Ho = 75 /175 km s~^ Mpc~^, and 
qo=0.5. 



2. SELECTION OF THE CLUSTERS 

Each target cluster in our sample was required to have 
bright (Lx > 5 x 10^^ h^^ ergs s~^), extended X-ray emis- 
sion in ROSAT observations. Therefore, the majority of 
the target galaxy clusters in our sample fall within the low 
redshift {z < 0.1) part of the X-ray-brightest Abell-type 
clusters of galaxies (XBACs, Ebeling et al. 1996) sample. 
The XBACs clusters are limited to high Galactic latitudes 
(|b| > 20°), low redshifts {z < 0.2), and ROSAT 0.1-2.4 
keV band X-ray fluxes above 5.0 x 10~^^ ergs cm~^ s~^. 
A further selection constraint was that each cluster was re- 
quired to have at least one linearly polarized radio source 
viewed through the X-ray emitting gas. Such radio sources 
are referred to as the cluster sample. A second set of po- 
larized (control) radio sources viewed, in projection, out- 
side the boundary of the X-ray emission was also selected 
for each target galaxy cluster. All polarized radio targets 

^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc 

^The ROSAT Data Archive is maintained by the Max-Planck-Institut fiir extraterrestrische Physik (MPE) at Garching, Germany. 



{cluster and control) were selected from the NRAO VLA 
Sky Survey (NVSS, Condon et al. 1998) data. 

Galaxy cluster selection was further constrained such 
that the sightlines of the target radio sources probed, col- 
lectively, the largest possible range of impact parameters. 
The sources also had to have sufficient polarized flux den- 
sity that follow-up polarimetry could be undertaken in 
short integrations (~ 5 minutes) at the VLA^. Speciflcally 
the inclusion criterion required Iia > 100 mJy and 1.4 
GHz polarization, mi. 4, greater than 1% in the NVSS sur- 
vey. A complete description of the galaxy cluster and radio 
source samples will be presented in Clarke, Kronberg, & 
Bohringer (in preparation). 

The final sample consists of 16 Abell clusters, 13 of 
which are members of the XBACs sample. The three 
non-XBACs clusters fall slightly below the flux limit for 
inclusion in the XBACs sample. Our 16 cluster sample 
was reduced from an original 24 that fell within the above 
criterion as severe radio frequency interference, which is 
endemic to some of the crucial (for RM) 20 cm bands, 
reduced the reliability of some RMs. To optimize data 
quality, we reduced the final sample to 16, since even this 
smaller number was statistically adequate. 

It is important to mention ab initio two types of sys- 
tematic bias that might result from the above selection 
criteria. First, the condition mi. 4 > 1% may preferen- 
tially select against regions of very high Faraday rotation, 
whose signature woiild be low polarization at the longer ra- 
dio wavelengths. This could have caused some high RMs, 
but not low RMs, to have been missed in our cluster sam- 
ple. This would statistically understate the clusters' true 
rotation measure distribution, and hence magnetic field 
strengths. Second, the very innermost regions of the clus- 
ter cores, which in some cases may have a cooling fiow, 
will have been missed because of their small angular cross 
sections. We do not consider this second form of bias to 
be serious, since this investigation is targeted to cluster 
volumes that do not have strong cooling flows. 

3. OBSERVATIONS AND ANALYSIS OF THE DATA 
3.1. Radio 

Target radio sources selected from the NVSS survey 
were re-observed with the VLA at four to six wavelengths 
within the 20 cm and 6 cm bands. These wavelengths were 
selected to provide Faraday rotation measures that are un- 
ambiguous within the range |RM| < 2600 rad m^^. The 
observations were undertaken in October and December 
1995, August 1996, and September 1997 in the VLA's B, 
D, and CS configurations respectively. 

The radio data were reduced within the NRAO AIPS 
package following the standard Fourier transform, decon- 
volve, and restore method. In addition, self-calibration 
was applied to each source to further reduce the effects of 
phase fluctuations. Images in the Stokes I, Q, and U pa- 
rameters were produced for each source at each of at least 
four wavelengths. 



3.2. X-ray 
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X-ray observations of each galaxy cluster were retrieved 
from the ROSAT Data Archive^. Thirteen of the target 
clusters were in the Pointed Observation archive while the 
data for the remaining three clusters were extracted from 
the ROSAT All-Sky Survey archive. All extracted X-ray 
data were taken using one of ROSAT's Position Sensitive 
Proportional Counters (PSPCs) which have moderate an- 
gular resolution (FWHM ~ 25") and are sensitive to pho- 
tons in range of 0.1 - 2.4 keV. 

The ROSAT X-ray data were reclucc^d using the Ex- 
tended X-ray Scientific Analysis Software (EXSAS, Zim- 
mermann et al. 1994) package under the European South- 
ern Observatory's (ESO's) Munich Image Data Analysis 
System (MIDAS). A radial X-ray surface brightness profile 
was determined for each cluster by integrating the PSPC 
photon events over concentric annuli of 15" and 30" for 
the ROSAT Pointed and RASS Observations respectively. 
The surface brightness profiles were fit with a hydrostatic 
isothermal model (Sarazin 1986). 

4. RESULTS 

The observed Faraday rotation measure of an individ- 
ual source is an algebraic sum of Faraday contributions 
due to our Galaxy, the cluster, any source-intrinsic com- 
ponent, and the general IGM. The latter three are usually 
small, and the Galactic RM contribution was statistically 
removed by subtracting the mean RM over all non-cluster 
sources within 10° of the cluster center from that of the 
radio probe in question. Due to the high Galactic latitudes 
of the target clusters, the mean Galactic contribution in 
the present sample is fairly small, on average 9.5 rad m~^. 
The cluster radio sources and associated Galaxy-corrected 
RMs are Usted in Table 1. 

In Figure 1 we plot the residual rotation measure of the 
radio sources as a function of cluster impact parameter 
in kiloparsecs. This figure displays a clear Faraday excess 
in radio sources viewed throTigh the X-ray emitting ICM 
(open points) as compared to those viewed beyond the 
detectable edge of the thermal cluster gas. The RM dis- 
tribution of the control sources has a width of 15 rad m~^ 
while that of the cluster sources viewed through the ICM is 
much broader at 114 rad m~^. The Kolmogorov-Smirnov 
test rejects the null hypothesis that the two samples were 
drawn from the same population with a confidence level 
of 99.5%. This confirms the detection of an intracluster 
Faraday rotating medium. 

Determination of magnetic field strengths at various im- 
pact parameters within our cluster sample requires some 
assumption about the field topology along the line of sight 
to the radio probe. The simplest model of the ICM mag- 
netic fields is a "uniform slab" in which the magnetic field 
has constant strength and direction through the entire 
cluster. Using this model and the X-ray determined elec- 
tron densities. Equation 1 yields magnetic field strengths 

1/2 

between ~ 0.5 3.0 /iG across the Faraday sample. 
More realistically, there are reversals along the line of sight 
to the radio source probe. In the simple case of an in- 
tracluster medium composed of J\f cells of uniform size 
and field strength, but random field directions, the field 
in an individual cell increases as vCV° over the uniform 
slab estimate. Using a simple tangled cell model with 
a constant coherence length, £ = Wh^^ kpc, yields an 



average intracluster magnetic field strength estimate of 

~ 5(^/10 kpc)~^/^ /Z75^ fiG. This coherence length £, esti- 
mated from RM images of three extended radio sources in 
our sample (sec below), is limited by the resolution of our 
images. Because < £(r) > may systematically change {e.g. 
increase with r) this "global" average field value could in- 

1 /2 

crease to values as high as 10 hj^ fiG near the cluster 
cores. 

The distribution of excess RM's appears to cutoff close 
to the observed X-ray outer boundary. This result, though 
interesting and new, requires more detailed X-ray and ra- 
dio data to understand the variation of the magnetic to 
thermal energy density esir) / eth{r) throughout the ICM. 

A striking effect, seen in Figure 1, is the virtual ex- 
clusion of small RMs at r < 500/if5^kpc. This sug- 
gests that the RM filling factor in normal galaxy clus- 
ters is very high. An independent, quantitative esti- 
mate of the RM filling factor is provided from analysis of 
multi-frequency polarization images of three extended ra- 
dio sources, 0039+212, 0056-013, and 1650-^815 (J2000) 
that are embedded within three of our clusters (Abell 
75, 119, and 2247 respectively). These sources project 
a combined area of 2.5 x lO^kpc^, and consist of 3 sets 
of 30, 48, and 24 contiguous independent RM sightlines, 
each of which has approximately the same areal cross- 
section: ^ (5kpc)^. The individual RM histograms across 
the three sources are consistent with normal distributions 
with means of -60, -144, and -97 radm~^. We find that 
95% of the sightlines within this subset of three clusters 
have RMs significantly higher than the dispersion for non- 
cluster sightlines, |RM| = 15radm~^. This implies that 
the areal filling factor of magnetic fields is at least 95% in 
the ICM. This strongly suggests that these enhanced mag- 
netic field levels permeate the clusters with a high filling 
factor, since within cell sizes down to a resolution of 10 
kpc almost no ray passing through the ICM escapes some 
magnetized region. 

5. DISCUSSION 

These results confirm the widespread existence of mag- 
netic fields in the central regions of non cooling flow clus- 
ters. The cluster-enhanced RM can generally be traced 
out to the periphery of the ROSAT-detcctable ICM X-ray 
emission. The rotation measure distribution across our 
cluster sample drops from ~ 200 rad m~^ in the central 
regions (which may be an underestimate, see §2) to the 
background level of ~ 15 rad m~^ at large radii. 

Our new measurements of (1) the RM, (2) the in- 
tracluster electron density, (3) the magnetic field vol- 
ume filling factor, and (4) the average tangling scale 
of the field enable us to estimate, even if only crudely, 
an important physical quantity — the total energy in 

the ICM magnetic field for "normal", non-merging, re- 

1/2 

laxed clusters. For a 5 /lyg /iC magnetic field in the 
inner 500/1^5^ kpc sphere, the total magnetic energy is 
Eb = 1.5 X 1061(^)3(^)2 ergs. The mag- 
netic energy content of the ICM can then be compared 
to the total thermal energy content in the same clus- 
ter volume. The latter is (again taking constant values 
within a fiducial radius that is close to both the RM 
cutoff radius in Figure 1 and the X-ray radius) Eth = 
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6.4 X Uf^n,, (^^ooT^) (to^-k) Kb ergs, where Ue is the 

intracluster electron density in units of 10~^ h^^^ cm"'^. 
The ratio is of order 2.5%, and the possibility that 
\B\ may be underestimated due to limited radio resolution 
makes it possible that could be even higher. Even 
at the lower bound of 2.5% the ratio suggests that the 
magnetic energy provides a non-negligible fraction of the 
energy budget of the ICM. 

We can now compare this approximate magnetic energy 
estimate with other sources of energy that are relevant for 
a cluster: the thermonuclear (stellar) energy released in all 
cluster-member galaxies, the; gravitational binding energy 
released from AGN/accretion disks, the cluster gas bind- 
ing energy, and the energy associated with past merger 
events. The available energy from stellar sources cannot 
be much greater than ~ 10*^^ ergs (Volk & Atoyan 2000), 
i.e. of comparable magnitude, and is thus insufficient to 
maintain the magnetic fields, barring a very high energy 
conversion efficiency. This means that thermonuclear en- 
ergy can be ruled out as the primary source of intracluster 
magnetic field energy. It must therefore be ultimately de- 
rived from gravity. A single powerful AGN/accretion disk 
can be expected to inject approximately 10^"'^ ergs over its 
lifetime into the ICM. Comparing the lifetime of the ra- 
dio source with the cluster lifetime we expect 10^ sources 
to have injected a total of ~ 10^^ — 10^'^ ergs into the 



ICM. This makes AGN/accretion disks an attractive pos- 
sible source of the ICM field energy, as has been suggested 
by Colgate & Li (1999). Other sources of ICM magnetic 
energy are the gravitational binding energy of the cluster 
gas, which is of order 10^^ ergs, and the energy associ- 
ated with a past cluster merger event, which is of order 
•|^q63-64 gj.gg This suggests that the magnetic energy pos- 
sibly "taps into" a combination of the gravitational energy 
released in AGN/accretion disks over the lifetimes of clus- 
ters, and shearing and shocks associated with larger scale 
infall of matter as the clusters evolve. 

Given that the cooling flows represent a late stage of 

1/2 

cluster evolution and our "global" value of 5 /iG for 
ICM zones within the clusters, but outside of cooling flow 
zones, we conclude that the cooling flow develops out of a 
medium whose field strength is already a significant frac- 
tion of what is seen in the cooling flow zones. 
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Table 1 

CLrsTER Sample Radio Sorru'Es 



Source 


Abell 




b b 


RM 


± RM 


(J2000) 




(lO'*'^ ergs s"i) 


(kpc) 


(rad m 




0039+212 


75 


0.26 


100 


-62.62 


9.03 


0040+212 


75 


0.26 


480 


-34.00 


10.95 


0042-092 


85 


3.96 


960 


4.88 


9.77 


0056-013 


119 


1.55 


280 


-149.16 


10.68 


0057-013 


119 


1.55 


1020 


14.02 


10.17 


0126-013a 


194 


0.06 


140 


94.22 


12.25 


0154+364 


262 


0.28 


570 


-202.36 


4.31 


0245+368 


376 


0.57 


630 


-48.06 


5.92 


0257+130 


399 


2.60 


280 


-185.89 


6.47 


0318+419 


426 


6.36 


550 


6.17 


15.40 


0316+412 


426 


6.36 


740 


74.94 


9.26 


0434-131 


496 


1.52 


310 


52.91 


6.69 


0434-133 


496 


1.52 


360 


35.89 


6.73 


0709+486 


569 


0.02 


3 


-229.74 


7.75 


0909-093 


754 


3.98 


1010 


-20.47 


9.16 


0908-100 


754 


3.98 


1460 


12.35 


9.43 


0919+334 


779 


0.07 


560 


16.55 


8.38 


1037-270 


1060 


0.23 


410 


9.38 


6.18 


1037-281 


1060 


0.23 


510 


25.21 


8.19 


1039-273 


1060 


0.23 


530 


103.99 


6.33 


1039-272 


1060 


0.23 


620 


103.73 


11.51 


1036-267 


1060 


0.23 


640 


25.76 


16.23 


1133+489a 


1314 


0.20 


270 


68.05 


6.42 


1133+490 


1314 


0.20 


310 


-50.38 


6.40 


1145+196 


1367 


0.68 


170 


257.46 


11.73 


1650+815a 


2247 


0.06 


200 


-127.36 


8.56 


1650+815b 


2247 


0.06 


270 


-131.75 


8.56 



"X-ray luminosity in the ROSAT 0.1 - 2.4 keV band determined from the current work. 
''Cluster-centric impact parameter of radio source. 
■^Galaxy-corrected rotation measure. 
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Fig. 1. — Galaxy-corrected rotation measure plotted as a function of source impact parameter in kiloparsecs for the sample of 16 Abell 
clusters. The open points represent the cluster sources viewed through the thermal cluster gas while the closed points arc the control sources 
at impact parameters beyond the cluster gas. Note the clear increase in the width of the RM distribution toward smaller impact parameter. 



